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Abstract 

This is essentially the text of my talk on Fiedler-Ptak scaling in max al- 
gebra delivered in the invited minisymposium in honor of Miroslav Fiedler 
at the 17th ILAS Conference in Braunschweig, Germany. 

Keywords: Max algebra, diagonal similarity scaling, nonnegative ma- 
trices, commuting matrices, matrix powers 

A MS codes: 15A80, 15B48, 15A27, 15A21 

0. Introduction 

The previous talks brought to our attention many ideas of Miroslav Fiedler 
that served as gold mines in mathematics. This talk will be based on only 
one gold mine: the Fiedler-Ptak diagonal similarity scaling and the role of this 
scaling in max algebra. 

We start from the initial appearence of this diagonal similarity scaling in a 
work of Fiedler-Ptak^ll , which was not meant to be applied in max algebra (an 
area which was only starting to emerge at the time of jTTj ) . 

We proceed with an introduction into max algebra and a specific version of 
the Fiedler-Ptak scaling appearing there, namely the max eigenvector scaling. 

Finally we concentrate on applications of the Fiedler-Ptak scaling in max al- 
gebra, in particular, to commuting matrices and matrix powers. In the following 
we abbreviate this to FP scaling. 

1. FP scaling Let Q(A) be a finite digraph weighted by ay S R+. Then, 
the FP scaling means existence of a potential y with positive coordinates yi such 
that y~ l a,ijyj < 1, which occurs if and only if \\ a,ij < 1 for any cycle C in 

Q{A). In the strong FP scaling the non-strict inequalities are replaced by the 
strict ones. 

Fiedler and Ptak used this scaling to characterize Hadamard preservers of 
certain classes of real- valued matrices. For a real- valued matrix A, define H(A) 

h u = \a u \, hij = -\dij\ for i ^ j. 
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Then class HI contains all matrices A such that H (A) has all principal minors 
positive, and class H», in addition, requires that the matrices have all non-zero 
diagonal entries. Here is one of the equivalence theorems (rather, a part of it) 
characterizing the Hadamard preservers of H and . 

Theorem 1 (Fiedler, Ptakp], Theorem 4.1) For B = (6 y ) E W nxn (where 
n>2), the following are equivalent: 

1. vfci, . . . , k s , i n- =1 b ^ +1 i < i nu h iki i + o; 

2. 3 diagonal positive D s.t. C = D~ 1 BD satisfying ^ \cu\ > \cij\ Vi, j G 
{l,...,n}; 

5. VAeH,Ao5e H; 

^. V4ei„ AoBeH,. 

FP scaling is used to establish equivalence between the first two statements. 

2. Max-algebra and FP scaling Consider nonnegative numbers R+ 
equipped with the ordinary multiplication ah and the new addition a © b := 
max(et, h). These operations are extended to matrices and vectors (of compatible 
sizes) in the usual way: we have {A®B)ij = (ay©6y), {A®B)ij = fe aifc<g>&fc.j 

k 

and A k = A ® • • ■ ® A for A = (a y ) and B = (6y). 

There are also max-plus and min-plus versions of max algebra. The first one 
is K U {— oo } equipped with a®b:=a + b and a b := max(a, b), the second 
one is M. U {+00} equipped with a (8) b := a + b and a © 6 := min(a, b), Both are 
isomorphic to our principal max-times version. 

Note that the max-algebraic matrix powers can be expressed in terms of 
optimal paths of fixed length. For a path P = i\ —>•••—>■ i}., denote the 
weight and the length of P by 

w(P) := a ili2 • . . . ■ a ik _ lik , l(P) = k - 1. 

Then the i,j entry of A* can be expressed as the greatest weight over all paths 
with length t, connecting i to j: 

of) = max{w(P) : i j, 1{P) = t}. 

Dropping the restriction on length leads to a max-algebraic analogue of (I— A)^ 1 
called Kleene star: 

A* := I © A © A 2 © • • • 

Thus a*j is the greatest weight of all paths connecting i to j (for i 7^ j). 
These obsevations also give a hint of applications of max-linearity in optimiza- 
tion and scheduling. 

The Kleene star converges if and only if there are no cycles C with w(C) > 1. 
In this case, A* = I © A © • • • © A" -1 . 
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In the following we denote by span ffi (A) the max-algebraic column span of 
A, i.e., the set of all max-algebraic combinations of the columns. In fact, all FP 
scalings can be described as positive vectors x in span ffi (A* ) , and all strong FP 
scalings as vectors x in the interior of span^ (A* ) . For a FP scaling, take a max- 
linear combination of all columns of A* with all positive coefficients, and for a 
strong FP scaling (possible only if span e (A*) is full-dimensional), take an ordi- 
nary linear combination. Actually, span^(j4*) is both a max-linear space (max 
cone) and an ordinary convex cone, and the interplay between these properties 
has been considered in [53], and Joswig-Kulas [IB]- See also Butkovic [5] where 
the strong FP scaling is used in the context of strong regularity of matrices in 
max algebra. 

Armed with this description of all FP scalings, Butkovic and Schneider [7] 
gave full solutions to several matrix scaling problems. One of these problems is 
to find all diagonal scalings X such that B = X~ x AX has equal column and 
row maxima: ba = maxj bij = maxj bji for all i. The solution to this problem 
is given by all positive vectors of span^(Q*), where Q := AD' 1 © D' 1 A and 
D = diag(an, . . . ,a„„). 

Another problem is, for A t ,B t ,C t £ R" x ™ such that g(A,) C Q(B % ) C g{d), 
where i = 1, . . . , m, to describe all X such that Ai < X~ 1 BiX < Ci. Informally, 
we have to get Bi sandwiched between Ai and Cj simultaneously by means of 
one diagonal similarity transform. To solve the problem, take positive vectors of 
span ffi (Q*), where Q 0™ x B.jC, ® 0™ 1 Af/Bf. Here, the divisions BJC, 
and Af /Bj are defined entry wise. 

Max eigenvector scaling is also a kind of FP scaling. Max eigenvectors are 
solutions to the spectral problem A® x — Ax. The greatest eigenvalue can be 
found as the maximum cycle geometric mean 



Further, the critical cycles are defined as the cycles where X(A) is attained, and 
the critical graph consists of all nodes and edges on the critical cycles. When 
X(A) = I (which can be assumed w.l.o.g.), the space of eigenvectors associated 
with X(A) can be explicitly described as the space (max cone) spanned by the 
columns of A* whose indices belong to the critical graph of A. 

Max eigenvector (with eigenvalue 1) corresponds to a scaling satisfying the 
following properties: 1) for all i,j we have x~ X dijXj < 1 , 2) for all i there is j 
such that x~ 1 aijXj — 1. (Also note that £ G C (A) implies x~ 1 a>ijXj = 1 for 
any FP scaling.) 

For any FP scaling x we can define the saturation graph Sat(A, x) consist- 
ing of the edges where x~ 1 a,ijXj = 1. We obtain that Q C {A) C Sat(A, x) 
and that the cycles of any saturation graph are precisely the critical cycles. 
Property 2) of the eigenvector scaling means that in the saturation graph of an 
eigenvector, each node has an outgoing edge. 

Below we give a schematic view of critical graph, and saturation graph of an 
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Figure 1: Critical graph and saturation graph 




Figure 2: The statement of Theorem [5] 



eigenvector scaling. 

3. Commuting matrices. When two matrices A and B commute, the 
eigenspace of one matrix A (corresponding to some eigenvalue) is invariant under 
the action of the other matrix B. Using this it can be proved that when A and 
B are irreducible, the intersection of eigenspaces always contains a positive 
common vector |17) . Take this vector x and use it for a simultaneous scaling of 
A and B. Further idea is to define A^ and ZjM as associated Boolean matrices 
of the saturation graphs Sa,t(A,x) and Sat(J3,x), and to observe that matrices 
A^ and £?M also commute. This leads us to study commuting digraphs 
with nonzero outdegree of each node. As an application, we may obtain 
properties of critical graphs of commuting max-algebraic matrices (since the 
strongly connected components of a critical graph are the same as the strongly 
connected components of any saturation graph). 

Commuting digraphs can be thought of as a "game" where two players (Alice 
A and Bob B) act together to connect one node to the other. If one node can 
be reached from another using some number of actions of A and B, then these 
actions can be arranged in every possible pattern, but the connecting routes 
(i.e., intermediate nodes and edges) will be different in general. Playing with 
such commuting digraphs we obtained the following 

Theorem 2 (Katz, Schneider, Sergeev [17]) Let Q\ and Qi commute and 
have nonzero outdegree of every node. Let Nh, ji = l,...,TOi, (resp. N„, 
v = 1, ...,m2 J be the node sets of nontrivial s.c.c. of Q\ (resp. of ^2)- Then 
there exists a cycle ci s Q\ with all nodes in (J N£ (resp. a cycle Ci £ G2 with 
all nodes in (J N*). 

Proof, (in pictures) Here is, essentially, what we need to prove: that there is 
a cycle of A (red) living in the s.c.c. "houses" of B (blue spheres). See Figured 
To show this, consider a subgraph of A living in a s.c.c. of B. If it has 
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Figure 3: The proof diagram 



cycles, then we are done. Otherwise, we consider a node where it has to go out, 
and we need to show that the destination set, i.e., the set of nodes to which 
A can move after it goes out, intersects with another s.c.c. of B. 

Now look at Figure [3] from left to the right (the destination set is the dashed 
ellipse). We ask B to turn around, since in the s.c.c. of B each node lies on a 
cycle, and then allow A to go out (left). But due to commutativity, we achieve 
the same effect if we first allow A to go out, and then ask B to connect to 
the destination node (middle). This implies that each node of the destination 
set has a predecessor in the graph of B, and hence there are cycles of B going 
through the destination set (right). Cycles of B belong to the strongly connected 
components of B, which finishes the proof. ■ 

4. Matrix powers. A theoretical motivation to study max-algebraic ma- 
trix powers is that they are similar to powers of graphs well-known in combina- 
torics [3]. In particular, the following fact has been known for a very long time, 
due to Cohen, Dubois, Quadrat et Viot [8]. 

Theorem 3 (Cyclicity Theorem) Let A e M" xn be irreducible nd \(A) = 
1, then the sequence A,A 2 ,A 3 ... is ultimately periodic. That is, after some 
coupling time T(A) there is a period 7 such that A r+1 — A r for all r > T(A). 

See also monographs on max algebra, such as El [15], for proof and more 
explanation, as well as [H QU E3 H3 El US for further development. 

Let us define the critical matrix A^- c ' as the associated matrix of the critical 
digraph. Then it satisfies (A k )^ = {A^) k . 

Assume that X(A) = 1 and that A is visualized (i.e., after an application 
of FP scaling). Let C, resp. R be the matrices extracted from the critical 
columns, resp. rows, of (A 1 )*, where 7 is the cyclicity of Q C (A), see Figure HI 
Let S := A^ (critical matrix). 

Then the following form of cyclicity theorem can be seen: 

Theorem 4 (Schneider, 2009) Let A E K" xn be irreducible and visualized, 
X(A) = 1, then there exists T(A) such that A 1 = CS l R for all t > T(A). 

In particular we obtain that A t+1 = A* for all t > T(A), S l R yields critical 
rows of A 1 and C«S" yields critical columns of A 1 for t > T(A). 

The matrix powers can be expressed in terms of optimal paths, and the 
cyclicity theorem also has a path sense. Namely, one defines strong paths as 
paths traversing a node of the critical graph. 
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Figure 4: Extracting C and R terms 



Theorem 5 (Sergeev, Schneider [23j ) Let Iif- r t be the set of strong paths 
connecting i to j of length t and let w(IIfj^) be the greatest weight of a path in 
this set. Then w(IL^ t ) = {CS t R) l:j after t> in 2 . 

We remark that the bound 3n 2 can be reduced to 2n 2 using the techniques 
of Hartman-Arguelles jT4j Theorem 4. 

In the scheduling applications it is important to bound T(A) from above. As 
a vague idea of doing this, let us start with a conjecture that after t > 0(n 2 ), 
we rather have a CSR expansion of the form 

A* = 0A*C fc S£i? fc , (1) 
k 

since at early times the matrix powers are not yet dominated by the critical 
graph and the main CSR term. 

An example of CSR expansion is shown below. This expansion is called the 
Nachtigall expansion, see QH [201 US • 

Namely, after formation of the main term we cancel all rows and columns 
with indices in the critical graph, and the second term is formed exactly in the 
same way. This procedure goes on until we get an acyclic matrix at some stage, 
and leads to a CSR expansion ([J). A bound on T(A) can be obtained using ([T]) 
and following the arguments of Hartman-Arguelles [M] Theorem 10: 

m , „ , max log an — min log a,,- 

T(A) < 2n 2 -5-^ — 6 13 . 2) 

log Ai - log A 2 

The Nachtigall expansion may seem too straightforward. A different scheme 
can be suggested using another kind of FP scaling called max-balancing. 
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Figure 5: Formation of first two terms in a Nachtigall expansion: a schematic 
example. 

Theorem 6 (H. Schneider, M. H. Schneider [21]) For any irreducible A G 
R" x " there is x G R" s.t. B = X AX has the following equivalent properties: 

• (max-balancing) For any cut M C {1 , . . . , n) , M := {f , . . . , n}\M , the 
greatest weight of the edges from M to M equals the greatest weight of the 
edges from M to M . 

• (cycle cover) For each edge there exists a cycle c containing this 
edge, where the weight of is minimal. 

An idea due to Hartman-Arguelles [14] is that in a max-balanced matrix, we 
can consider a threshold digraph consisting of all edges with weights not pro- 
ceeding some parameter. Lowering this parameter until the threshold digraph 
has a new strongly connected component determines the "region of influence" 
of the main CSR term, and in the same vein, of all the rest to be formed. This 
idea will be made precise in a forthcoming work of the author. 

5. Other works. 

The author's works on application of FP scaling were inspired by the study of 
power method in max algebra due to Eisner and van den Driessche [5] . We also 
mention that the same authors use FP scaling in their study of symmetrically 
reciprocal matrices [10] . 

The idea of FP scaling was further extended to semigroups by Gaubert [12] 
who gave a positive solution to Burnside problem for the semigroups of matrices 
in max algebra. In a more recent work, Benek Gursoy and Mason [3] extend the 
strong FP scaling to the set of matrices, applying it to study the asymptotic 
stability of max-algebraic difference inclusions. 

The same idea was crucial for the development of max-algebraic methods in 
the (ordinary) eigenvalue perturbation theory, aimed to improve and extend the 
theory of Visik, Ljusternik and Lidskh. Namely, Akian, Bapat and Gaubert [I] 
consider matrices of the form 

(A e )ij = a l0 e A ^ + o{e A *>), e ^+ 

where a G C nx ™, A G (RU {+oo})™ xn . One of the goals is to give the first order 
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asymptotics of the eigenvalues 



(3) 



It can be observed that the ordinary FP scaling (of the min-plus matrix A, in 
terms of min-plus algebra) already helps, visualizing some entries of a (corre- 
sponding to the critical graph of A) which are really important for the asymp- 
totics. Further idea is to extend the critical graph to cover the whole A, which is 
achieved by means of a diagonal similarity scaling similar to max-balancing, to 
leave out all unimportant entries of a and to construct a sequence of Schur com- 
plements in what remains of a. This sequence is determined by a critical graph 
tower associated with A (in min-plus algebra). This critical graph tower also 
determines the exponents Aj in ([3]), and appear as groups of eigenvalues (in 
the ordinary linear algebra), determined by the sequence of Schur complements. 
See Q] Theorem 5.1 for precise formulation. 

We also mention that the diagonal similarity scaling of FP type is useful 
in the idempotent analysis (extension to max-algebraic functional spaces), see 
Kolokoltsov and Maslov [TH] Chapter 2. 
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